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Abstract: A method is proposed to estimate the velocity field of an unsteady flow using a limited 
number of flow measurements. The method is based on a non-linear low-dimensional model of the 
flow and on expanding the velocity field in terms of empirical basis functions. The main idea is to 
impose that the coefficients of the modal expansion of the velocity field give the best approximation 
to the available measurements and that at the same time they satisfy as close as possible the non- 
linear low-order model. The practical use may range from feedback flow control to monitoring of 
the flow in non-accessible regions. The proposed technique is applied to the flow around a confined 
square cylinder, both in two- and three-dimensional laminar flow regimes. Comparisons are provided 
with existing linear and non-linear estimation techniques. 
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Un observateur non lineaire pour des ecoulements 
tridimensionnels instables 



Resume : Nous proposons une methode pour estimer le champ de vitesse d'un ecoulement instable 
en utilisant un nombre limite de mesures. La methode est basee sur un modele d'ordre reduit non 
lineaire de l'ecoulement et sur l'expansion du champ de vitesse en termes de fonctions de base em- 
piriques. L'idee principale est d'imposer que les coefficients de l'expansion modale du champ de 
vitesse donnent la meilleure approximation des mesures disponibles et qu'en meme temps ils satis- 
fassent aussi etroitement que possible le modele d'ordre reduit non lineaire. L' utilisation pratique 
peut s'etendre du controle actif en boucle fermee de l'ecoulement a la surveillance de l'ecoulement 
dans des regions inaccessibles. La technique proposee est appliquee a l'ecoulement autour d'un cy- 
lindre carre confine, dans des regimes laminaires bi et tridimensionnels. Des comparaisons avec des 
techniques lineaires et non lineaires d'evaluation existantes sont fournies. 

Mots-cles : modeles reduits, estimation dynamique, observateurs non lineaires 
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1 Introduction 

The problem of deriving an accurate estimation of the velocity field in an unsteady complex flow, 
starting from a limited number of measurements, is of great importance in many engineering ap- 
plications. For instance, in the design of a feedback control, a knowledge of the velocity field is a 
fundamental element in deciding the appropriate actuator reaction to different flow conditions. In 
other applications it may be necessary or advisable to monitor the flow conditions in regions of space 
which are difficult to access or where probes cannot be fitted without causing interference problems. 

The method that we propose exploits an idea which is similar to that at the basis of the Kalman 
filter (see ifTTI ). The starting point is a Galerkin representation of the velocity field u(x, t) in terms 
of N r empirical eigenfunctions, «&*(a:), obtained by Proper Orthogonal Decomposition (POD) (see 



where u{x, t) : K" X [0, T] -» R", &(x) : R n -> M", n e {2, 3} according to the physical space 
dimension, u(x) is some reference velocity field and a;(i) : I = [0, T] Cl->1. 

For a given flow, the POD modes can be computed once for all based on Direct Numerical 
Simulation (DNS) or on highly resolved experimental velocity fields, such as those obtained by 
particle image velocimetry. An instantaneous velocity field can thus be reconstructed by estimating 
the coefficients (t) of its Galerkin representation. 

One simple approach to estimating the POD coefficients is to approximate the flow measure- 
ments in a least square sense, as done, for instance, in (9j- 

A similar procedure is also used in the estimation based on gappy POD, see lfl8l and lfl9l . 
Another possible approach, the linear stochastic estimation (LSE), is based on the assumption that a 
linear correlation exists between the flow measurements and the value of the POD modal coefficients 
(see, for instance, 13). 

However, these approaches encounter difficulties in giving accurate estimations when three- 
dimensional flows with complicated unsteady patterns are considered, or when a very limited number 
of sensors is available. Under these conditions, for instance, the least squares approach cited above 
(LSQ) rapidly becomes ill conditioned. This simply reflects the fact that more and more different 
flow configurations correspond to the same set of measurements. To circumvent those problems, 
many contributions in the literature have aimed to determine the effective placement of the sensors 
(see e.g. fl4l . (6), 0, |fl9l ). For example in |fl9l , a systematic approach to sensor placement is 
formulated within the gappy POD framework using a condition number criterion. 

In order to improve estimation performance, extensions of the above methods have been pro- 
posed: quadratic stochastic estimation (QSE) |fl~), fi"3ll and spectral linear stochastic estimation 
(SLSE) 1 7 1. They allow more accurate estimations compared with LSQ or LSE methods, but, in fact, 
neither of these methods takes into account the underlying dynamic model that the POD coefficients 
must satisfy, i.e., a finite dimensional equivalent of the Navier-Stokes equations that is obtained by 
the Galerkin projection of the flow equations on the POD modes retained for the representation of 
the velocity field. In this sense, the aim of the present study is to discuss an approach that combines 
a linear estimation of the coefficients di(t) with an appropriate non-linear low-dimensional flow 
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Figure 1: Computational domain f2. 



model. Our objective is not, however, to propose an estimation method that can be readily imple- 
mented for real time applications, even if a few indications in this direction are given. Rather, our 
objective is to understand whether a non-linear observer outperforms existing linear flow observers, 
without the constraints imposed by an actual recursive algorithm, e.g., a real-time computation. 
Moreover, instead of what was done, for example, in iflOl , this study is confined to a deterministic 
framework, since the model as well as the measurements are supposedly not affected by noise. If, 
within this framework, a dynamic estimation turns out to be less satisfactory than static approaches, 
i.e., those which use no model, then there would be little interest in pursuing the research in this 
direction. In addition, we address the issue of the sensitivity of the proposed approach to sensor type 
and location. Finally, we present an application to a flow, which is characterized by a significant 
three-dimensionality and a non-periodic dynamics. 

2 Flow set up and low order model 

The flow over an infinitely long square cylinder symmetrically confined by two parallel planes is 
considered. A sketch showing the geometry, the frame of reference and the adopted notation is 
plotted in figure Q] At the inlet, the incoming flow is assumed to have a Poiseuille profile with 
maximum center-line velocity U c . Two Reynolds numbers Re = U c L/v were considered, one at 
which the flow is two-dimensional (Re = 150) and the other one leading to a three-dimensional 
flow in the wake (Re = 300). With reference to figure Q] L/H = 1/8, Lin/L = 12, Lout/L = 
20 and Lz/L = 0.6 for the two-dimensional case, whereas L/Lz = 6 for the three-dimensional 
one. Periodic boundary conditions are imposed in the span-wise direction and no-slip conditions 
are enforced both on the cylinder and on the parallel walls. Details concerning the grids and the 
numerical set up are reported in 0. All the quantities reported in the following have been made 
non-dimensional by L and U c . The two-dimensional flow obtained at Re = 150 is a classic vortex 
street with a well defined shedding frequency. However, the interaction with the confining walls 
adds to the complexity of the flow and leads to some peculiar features, like the fact that the vertical 
position of the span-wise vortices is opposite to the one in the classic von Karman street (JH). For 
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the three-dimensional case the situation is even more complex, due to instabilities developing in the 
span-wise direction. The flow is no longer periodic and exhibits complicated flow patterns (15]). 
The POD modes $ k (x) are found using the snapshot method (03]) 



where — u(x, U) are flow snapshots taken at time ti G [0, T], N is the number of snapshots, 
k G {1, . . . ,N}, and the coefficients if e 1 are such that the vectors {b\, . . . , 6jy) are the eigen- 
vectors of the time correlation matrix J n U^' ■ U"' dx, of size N x N. Only a limited number of 
modes, N r , is used to represent the velocity field. In particular we took N r = 6 and N r = 20 for the 
two-dimensional and the three-dimensional cases, respectively. 

A Galerkin projection of the incompressible Navier-Stokes equations over the retained POD 
modes has been carried out. This leads to the following N r -dimensional dynamical system 



where a(t) : I — > M. Nr and a(t) = {ai(t), . . . , a,N r (t)}\ r, k and s run from 1 to N r and the 
Einstein summation convention is used. The scalar coefficients B ksr come directly from the Galerkin 
projection of the non-linear terms in the Navier-Stokes equations, and they can easily be expressed 
in terms of the POD modes. The scalar terms A r and C kr are calibrated using a pseudo-spectral 
method to take into account the pressure drop, as well as the interaction of unresolved modes in 
the POD expansion. The calibration consists in solving an inverse problem, where the coefficients 
A r and C kr are found in order to minimize the difference, measured in the L 2 norm, between the 
model prediction and the actual reference solution. See J8] for a detailed discussion of the calibration 
technique. 

The resulting model for the two-dimensional flow configuration considered here is very accurate 
in describing the asymptotic attractor ((9] and O). For the three-dimensional case, it was shown 
in UJ that the calibrated model is capable of accurately reproducing the complicated flow dynamics 
resulting from the interaction of the three-dimensional vortex wake with the confining walls inside 
the calibration interval. This is shown in figure |2j where the predictions of some POD modal co- 
efficients given by the dynamic POD model within the calibration interval are compared to those 
obtained from the projection of the fully resolved Navier-Stokes simulations. However, these results 
rapidly deteriorate as soon as the model becomes an actual prediction tool, i.e., outside the time 
interval in which the calibration is performed. Nonetheless, the model can be used in conjunction 
with experimental measurements in order to reconstruct the entire flow field, as will be shown in the 
next section. 

3 Non-linear observer 

Our aim is to provide an estimation of the modal coefficients Oi (t) starting from N s flow measure- 
ments fk, k € {1, . . . , N s }. Let a.i(t) be the projection of the velocity field u(t) over the i-th POD 
mode and oti it) be its estimated value at time t. 



R r (a(t)) = a r (t) -A r - C kr a k {t) + B ksr a k {t)a s {t) = 
a r (0) = (u(aj,0) -«(a:),* r ) 



(2) 
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Figure 2: Three-dimensional flow: projection of the fully resolved Navier-Stokes simulations over 
the POD modes (continuous line) vs. the integration of the dynamical system inside the calibration 
interval, obtained retaining the first 20 POD modes (circles). Only 8 representative modal coeffi- 
cients are shown here. 



We assume that each measurement f k is a scalar quantity which depends linearly on the instanta- 
neous velocity field u(t). For instance, fk can be a point- wise measurement of a velocity component, 
or of a shear stress, or it can be a spatial average of a linear combination of velocity components. 

The available spatial information may be exploited by using a LSQ approach, as done in |9). At 
any given time t, thanks to the linearity of f k with respect to u and to the modal decomposition of 
the velocity field (see Eq. (QJ), fk can be written in terms of POD modes 

N r 

f k ( U (r))c^2 aj (r)f k {^) (3) 

3=1 

where fk is obtained from the application of f k to the vector field associated to mode 
Then, the following least-squares problem has to be solved for every r 

N s I N r 
{oi(r),...,o JVr (T)}^ I ^ 

This problem leads to the solution of a N r -dimensional linear system of equations. Once this prob- 
lem is solved, the estimated modal coefficients can be written 

a j (r) = Y^T kj f k (u(r)) (5) 
fe=i 
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where T is a known rectangular matrix of size N s x N r , The error minimization (0]i leads to a linear 
representation of the estimated modes as a function of the measurements. 

The LSE approach, conversely, exploits temporal rather than spatial information and is based on 
the assumption that a linear relation exists between the modal coefficients and the measurements 

a j (T)=Y i A kj f k (u(T)) (6) 

k=l 

where A is now an unknown rectangular matrix of size N s x N r . This matrix is determined by 
imposing the condition that Mj G {1, . . . , N r } and Vfc G {1, . . . , N s } 

/ a 3 {t) f k (u (t))dt= V k mj f m (u (t)) f k {u (t)) dt (7) 
Jo Jo m=1 

The time interval [0, T] is the same as that considered for building the POD modes. Hence, since the 
left-hand side is known, a set of linear equations is obtained; these uniquely define the matrix A. 

The LSQ and and LSE both provide linear estimation of the modal coefficients. Matrices T 
and A have the same size, although the coefficients are different. In the following we overcome the 
assumption of a linear relation. 

Let us assume that a certain number of measurements at consecutive times r m , m G {1, . . . , N m } 
are available. The main idea of the dynamic-estimation approach proposed here is to impose that the 
coefficients of the modal expansion of the velocity field give the best approximation to the available 
measurements, using either LSQ © or LSE (|6|, and that at the same time they satisfy as closely as 
possible the non-linear low-order model (|2j. 

In the LSQ case this is done by minimizing the sum of the residuals of (0 and the residuals of 
(0 for all times r m . More precisely, let a(t) : K — > WL Nr and a(t) = {ai(t), . . . , ajv r (t)}, we have 

a(t) = argmin ^ C ' R ^T R 2 r (a{T m )) + ^(a r (r m ) 

a (*) m=l \ r=l r=l 

where a(t) = {ai(t), . . . , a,N r (t)}- The parameter Cr weights more the measurements (LSQ) or 
the dynamic model in the definition of the residual norm. It could be systematically tuned, or it could 
be a matrix. In the numerical experiments reported here, this parameter has been in a heuristic way, 
leaving a consistent analysis to future investigations. The minimization of this functional is reduced 
to a non-linear algebraic problem. As in J8j, a pseudo-spectral approach is used and each a r (t) 
is expanded in time using Lagrange polynomials defined on Chebyshev-Gauss-Lobatto collocation 
points. The necessary conditions for the minimum result in a non-linear set of algebraic equations for 
the coefficients of the Lagrange polynomials. The solution is obtained by a Newton method, which, 
in the present applications, usually converges in a few (typically 5 to 8) iterations. The solution of the 
problem ([8]l provides an estimation for the POD modal coefficients for all modes and for all instants 
at which measurements are available. This allows the reconstruction of the entire flow field at the 
same instants through equation (Q~|). The above method, therefore, represents a non-linear observer 
of the flow state. In the following, it will be referenced as K-LSQ. 



(8) 



fe=i 
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A similar approach can be obtained for the LSE technique, by substituting in Eq. (© the residuals 
of Eq. (O instead of those of Eq. (0. This approach is referenced as K-LSE. 

In literature, there exist other flow estimation techniques that are non-linear in the flow measure- 
ments. In the following we will compare the results of the proposed non-linear dynamic estimation 
to one of them, a quadratic extension to LSE (|fl~), fPUl ). This method is based on the assumption 
that equation (|6]l is just the first term of a Taylor expansion with respect to the sensor measurements, 
whereas QSE takes into account the second order term, too. Hence, we have 

JV S N, N s 

aj (f) = £ c« (t)) + 2 E ( M ( r )) fm ( M (9) 

fe=l fe=l m=l 

where the scalar coefficients A^j and tlkmj are obtained using double, triple and quadruple cor- 
relations between measurements in an equation equivalent to (0. This approach is referred to as 
QSE. 

Once the matrices appearing in equations ((5]) (O and (0 are computed, the estimation of the 
modal coefficient at a certain time is based on the measurements made at the same time. In contrast, 
10, lfl7l proposed to take into account integrated temporal data by assuming a linear dependence 
between the modal coefficients and the flow measurements in a non-local way, by working in the 
frequency domain. Let a be the Fourier transform of a and fj that of fj, then for each frequency we 
pose 

N, 

aj = Y^r kj f k do) 

fe=i 

where is a matrix obtained by appropriate ensemble averages and depends on the frequency. 
In the time domain this amounts to a convolution integral between the measurements and the time 
dependent matrix T. We call this approach SLSE. As compared to QSE and SLSE, the dynamic 
estimation procedure that we propose is non-linear and, at the same time, it takes into account the 
evolution of the modal coefficients in time by constraining such evolution to a model, in the weak 
sense determined by [8] 

Concerning the applicability of the methods described above, it is important to recall that the LSE 
and LSQ approaches are readily applicable to real-time estimation, as well as the QSE, although the 
cost of this last approach scales as instead of linearly as in the previous two cases. Conversely, 
the SLSE approach is more difficult to be used for real-time estimations, since it uses the whole 
temporal history of the measurements, collected in a time interval, coupled together (linearly) via 
the Discrete Fourier Transform (DFT). This implies that the estimation problem must be tackled 
after having collected enough temporal information and it consists of a number of LSE problems 
equal to the number of retained frequencies, plus additional DFT's of the measurements and of the 
estimated POD coefficients. 

Similarly to what done in the SLSE approach, in the present dynamic estimations the temporal 
histories of the measurements are coupled together (non-linearly) by the dynamic POD model. This 
aspect poses difficulties in a real-time application. Indeed, as pointed out in the introduction, K- 
LSQ or K-LSE are thought to be applied a-posteriori, because their computational cost, although 
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unimportant in a post-processing phase, is large for a real-time analysis. Nevertheless, although 
actual real-time applications are premature, a proposal for their prospective implementation for real- 
time estimation is the following. The flow state at a given time t* could be estimated by considering 
the measurements taken at that time and at the previous N m — 1 ones. At the successive sampling 
time, the corresponding new measurements are added and the oldest ones are dropped, keeping the 
number of measurements considered constant. In other words, reconstruction is carried out using a 
fixed number of measurements distributed in a time interval which is located before t*, and which 
translates as time increases. The sampling rate (i.e. r m — T m _i) and N m can be tuned in order 
to decrease the computational costs while granting the level of accuracy required by the particular 
application. Moreover, when a new set of measurements is added, the Newton method for solving 
the non-linear system would be restarted from the previous solution, which is already close to the 
final solution, thus definitely reducing the number of iterations for convergence. 

In contrast with the other methods, the proposed approaches need a working Galerkin model as 
a fundamental ingredient. The construction of such a model can be carried out from the information 
needed to build the POD database, a necessary step for all the methods considered here. Therefore 
no additional information is needed as compared to other approaches. 

4 Results and discussion 

The K-LSQ and K-LSE are used to reconstruct the flow in the configuration described in Sec. [2] 
both in the two- (Re = 150) and three-dimensional (Re = 300) cases. Results are compared to those 
obtained by the most common techniques available in the literature and cited in the Introduction. 

Accuracy in the prediction of the single modal coefficients and in the reconstruction of the veloc- 
ity fields were appraised. In both cases, differences with respect to the reference case (DNS) were 
quantified in terms of relative error in the L 2 norm, i.e., the L 2 norm of the difference between the 
estimated and the reference quantity divided by the norm of the reference quantity. 

Several parameters are involved in the set-up of the K-LSQ and K-LSE models. They are re- 
lated to (i) the dynamic POD model: number of retained modes, calibration interval, number and 
temporal distribution of available snapshots; (ii) the selected flow measurements: number, type and 
collocation. 

The parameters involved in the derivation of the POD dynamical system used herein are the same 
that were selected in previous studies (H, Q). As for the flow measurements, both velocity and 
shear-stress sensors were used. While velocity measurements are often considered in the literature, 
due to their widespread use in practice, shear-stress sensors are less common. Nevertheless, they 
were used here mainly because they are challenging from a numerical point of view, as they involve 
spatial derivatives of the POD modes. Also, they can be implemented in practice although limitations 
of accuracy and time resolution may exist (see, for instance, Ifl6l ). Different sensor locations were 
tested, to account for the sensitivity of the proposed approaches to sensor placement. Since the 
performance of the standard techniques such as LSE or LSQ is influenced by sensor placement, 
some sensor configurations were selected following the suggestions given for LSE in |6|. On the 
other hand, none of the considered sensor configurations is optimized for K-LSE or K-LSQ, in order 
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to verify the sensitivity of such methods with respect to sensor placement. In fact, optimal sensor 
placement may turn out to be a time-consuming operation for complex three-dimensional flows. 

4.1 Two-dimensional case: Re = 150 

For this rather simple flow, we consider the situation in which a limited number of measurements are 
available, i.e. only 2 sensors. Three different configurations were analyzed, two involving stream- 
wise velocity sensors and one involving shear-stress sensors. 

The velocity sensors were placed in relation to the spatial structure of the streamwise component 
of the first two POD modes. In particular, in the first configuration one streamwise velocity sensor is 
placed on the maximum of the first POD mode which is closest to the cylinder (PI ~ (2.39, 0.52)) 
and one in the middle between PI and the minimum of the second POD mode closest to the cylinder 
(~ (1.96, 0.50)). The second configuration has the first streamwise velocity sensor in PI and the 
second one in the point (1.98, —0.76). A third configuration was considered with two shear-stress 
sensors located on the confining walls (y = ±4.0) at x = 4, in a region which satisfies the following 
criteria on a shedding cycle: the rms value of the shear-stresses is maximum and the reconstruction 
error of the shear-stresses is minimum for a given number of POD modes. 

The POD low-order model of the two-dimensional flow is obtained by using 95 snapshots, uni- 
formly distributed throughout two vortex shedding cycles (T ~ 13 is the non-dimensional duration 
of the time interval), and by retaining iV r = 6 modes. The calibration of the model is performed 
in the same interval using 81 collocation points. As shown in (9), the calibrated model accurately 
reproduces the flow inside and outside the calibration interval. 

Figure [3] shows the POD modal coefficients predicted by K-LSQ and by LSQ considering the 
third (shear-stress sensors) configuration, together with their reference values found from DNS. In 
addition, errors in the prediction of the modal coefficients given by LSQ, LSE, QSE, K-LSQ and 
K-LSE in the first (velocity sensors) and third (shear-stress sensors) configurations are reported in 
Table [TJa) and (b), respectively. The values obtained for the second considered sensor configuration 
are not shown since they are very similar to those of the first one. 

The time interval over which reconstruction is performed is approximately 13 time-units long 
(non-dimensional time); it contains two shedding cycles, and it starts just after the end of the time 
interval on which the POD model was built and calibrated. 

Tables \V[a) and (b) show the relative reconstruction errors on the velocity components and on 
their fluctuating part. It appears that two (velocity or shear-stress) sensors are not sufficient to obtain 
reliable predictions of the modal coefficients by LSQ, LSE, even if LSE leads to a better estimation 
than LSQ. Accuracy problems persist also with the QSE approach, even if in this case the predic- 
tions are more accurate than those obtained with LSE. The results obtained by LSE and QSE are 
more accurate in the configuration with velocity sensors because their positions were selected to be 
effective for LSE, as already pointed out. Nevertheless, the errors on the estimation of the first two 
modal coefficients are still large. This leads to severe errors in the estimation of the fluctuating part 
of the velocity field since the first two POD modes represent about 94.8% of the fluctuating energy. 
Even if the mean flow energy is important with respect to the fluctuating energy, errors in the modal 
coefficients lead to detectable errors in the reconstruction of the velocity components. Note that the 
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Figure 3: (a) and (b): POD modal coefficients ai and estimated with the K-LSQ and the LSQ 
approaches, compared to their exact values computed by projection of the DNS velocity fields, (c) 
coefficient ae, DNS and LSQ; (d) coefficient a^, DNS and K-LSQ. Note that in (c) and (d) different 
axis scales are used. 

reconstruction errors on the vertical component are larger than those on the streamwise one simply 
because the contribution of the mean flow on that component is much lower than in the streamwise 
component. Table Q] shows that both K-LSQ and K-LSE give an accurate estimation not only of the 
first two modal coefficients, but of all the retained modes. This leads to a precise estimation of the 
velocity field as well as of its fluctuating part. Moreover, the accuracy of the results is very similar 
with shear-stress and velocity sensors, indicating a weak sensitivity of the approach with respect to 
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Tab. 


e(ai)% 


e(a 2 )% 


e(a 3 )% 


e(a 4 )% 


e(a 5 )% 


e(a 6 )% 


LSQ 


108.41 


340.48 


538.59 


2210 


10900 


8340 


LSE 


71.06 


27.12 


99.71 


97.97 


99.91 


99.93 


K-LSQ 


0.47 


0.55 


2.58 


2.66 


4.65 


4.67 


K-LSE 


0.82 


0.76 


9.82 


9.82 


14.98 


15.59 


Tab. Htb) 


e(ai)% 


e(a 2 )% 


e(a 3 )% 


e(a 4 )% 


e(a 5 )% 


e(a 6 )% 


LSQ 


100.33 


140.60 


752.10 


1040 


5490 


37800 


LSE 


46.87 


91.97 


102.83 


100.99 


101.40 


100.21 


K-LSQ 


0.06 


0.09 


6.06 


6.09 


9.72 


9.56 


K-LSE 


2.99 


3.08 


7.07 


8.31 


17.99 


18.47 



Table 1 : Relative percentage errors (in L 2 norm) on the estimation of the POD modal coefficients 
(e(cii)) in the first (a) and third (b) sensor configuration. In this case time-averaging is carried out 
over the estimation time period. 



Tab. Ha) 


e{U)% 


e{V)% 


e(U')% 


e(V')% 


LSQ 


106.79 


937.42 


906.96 


1370 


LSE 


6.32 


37.14 


53.96 


54.26 


K-LSQ 


0.63 


3.97 


5.39 


5.80 


K-LSE 


0.69 


4.42 


5.93 


6.46 


Tab.[2}b) 


e{U)% 


e(V)% 


e(U')% 


e(V')% 


LSQ 


273.22 


2490 


2350 


3640 


LSE 


8.05 


46.14 


68.32 


67.44 


K-LSQ 


0.65 


4.10 


5.54 


6.00 


K-LSE 


0.77 


4.91 


6.54 


7.17 



Table 2: Relative percentage errors (in L 2 norm) on the estimation of the velocity components 
(e(U),e(V)) and of their fluctuating part (e(U'),e(V )), in the first (a) and third (b) sensor configu- 
ration. In this case time-averaging is carried out over the estimation time period. 

the type and location of the sensors. This is not the case for the LSQ, LSE and QSE methods, which 
show a higher sensitivity to this aspect, confirming what has already been reported in the literature. 

We compare our results to those of |6), Tab. 2(a) 13-th case. With LSE and 2 sensors they 
found e(ai) ~ 76.6% and e{ci2) — 15.1%, errors that are similar to those reported in Table[TJa) for 
LSE. Using the dynamic estimation, the errors on the same coefficients are two orders of magnitude 
lower. Furthermore, using K-LSQ method and two shear-stress sensors (Table [Ub)) the first two 
modal coefficients are estimated with an error lower than 0.1%, i.e., three orders of magnitude lower 
than LSE. 
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4.2 Three-dimensional case: Re = 300 

The flow patterns in this case are definitely more complex than those of the previous one. For 
this reason, 24 flow measurements were used for the reconstruction procedure, organized in five 
different configurations, two involving only velocity measurements and three involving both veloc- 
ity and shear-stress measurements. In the last three configurations, the shear-stress sensors were 
selected following the same criterion adopted in the 2D case, they are 14 in number and in all 
considered cases are symmetrically placed on both the confining walls (y = ±4) at x = 4 and 
z = {1.2, 1.5, 2.7, 3, 3.3, 4.5, 4.8}. The placement of the velocity sensors has again been chosen on 
the basis of the spatial structure of the streamwise velocity of the first 12 POD modes. The different 
configurations are listed below, together with a brief description of the rationale for the placement 
of the velocity sensors: 

1 . 24 velocity sensors distributed on 6 equispaced slices in the axial (z) direction; on each slice, 
the sensors are on the lines connecting the maximum and minimum closest to the cylinder of 
the first two POD modes. On each segment, the sensors are approximately in the middle, but 
slightly closer to the extrema of the first POD mode. 

2. 24 velocity sensors distributed on 4 equispaced slices in the axial (z) direction; on each slice, 
3 points are selected in the region of overlapping between the maxima and minima of the 
low-frequency POD modes (modes 3, 4, 7, 8, 9 and 10) and 3 on the overlapping region of the 
extrema of the vortex shedding modes (modes 1, 2, 5, 6, 7, 11 and 12) (see f3| for details on 
the separation between low-frequency and vortex-shedding POD modes). 

3. 14 shear-stress sensors and 10 velocity sensors distributed on 5 equispaced slices in the axial 
(z) direction; on each slice, the velocity sensors are placed on the maximum and minimum 
closest to the cylinder of the first POD mode. 

4. 14 shear-stress sensors and 10 velocity sensors. 6 equispaced slices in the axial (z) direction 
are considered. On 4 slices, 2 velocity sensors are placed as in the previous case. Two sensors 
are placed on the remaining slices, corresponding respectively to the maximum and minimum 
of the third POD mode (low frequency mode). 

5. 14 shear-stress sensors and 10 velocity sensors located in the wake, at the points reported in 
table [3] 

The low-order model of the developed three-dimensional flow was derived retaining the first 20 POD 
modes obtained from a database of 151 snapshots, uniformly distributed over eight vortex shedding 
cycles (~ 52 non-dimensional time units). Calibration was carried out in the same time interval, 
and the results obtained integrating the dynamic model within the calibration interval are reported 
in figure[2] More details of this model can be found in [3]] (the model is denoted there as POD2). In 
contrast to the two-dimensional case, the POD model is inaccurate outside the calibration interval, 
so it cannot be used as a predictive tool by itself. The estimation of the flow is carried out here in two 
different time intervals, one starting just after the end of the time interval in which the POD model 
was calibrated and the other one being about 52 time units distant from the calibration interval. Both 
the intervals are approximately 30 time units long, including approximately 4 shedding cycles. 
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Velocity sensor 


X 


y 


z 


Velocity sensor 


X 


y 


z 


1 


5.01 


1.03 


2.00 


6 


5.01 


-1.03 


4.00 


2 


6.99 


1.04 


2.00 


7 


6.99 


-1.04 


4.00 


3 


6.03 


-1.04 


2.00 


8 


6.03 


1.04 


4.00 


4 


6.01 


-0.99 


2.00 


9 


6.01 


0.99 


4.00 


5 


5.96 


-0.97 


2.00 


10 


5.96 


0.97 


3.99 



Table 3: Positions of the velocity sensors in the three-dimensional case 



415 420 425 430 435 440 

Non-dimensional time 

Figure 4: Relative percent error in the reconstruction of the U component projected on the retained 
POD modes as a function of time, when varying Cr. 
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Figure 5: Estimation of some representative modal coefficients in the three-dimensional case, for 
the configuration (b) (see text) and in the time interval close to the calibration one, together with the 
reference ones evaluated from the DNS simulation. Note that different axis scales are used on the 
left and on the right plots. 
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SC (a) 


e(U)% 


e(V)% 


e(W)% 




SC (b) 


e(U)% 


e(V)% 


e(W)% 


KLSQ 
KLSE 


9.85 
9.60 


37.00 
36.49 


108.12 
103.92 


KLSQ 
KLSE 


10.00 
10.03 


36.65 
38.19 


111.52 
105.57 




e(U')% 


e(V')% 


e(W')% 




e(U')% 


e(V')% 


e(W')% 


KLSQ 
KLSE 


61.89 
9.60 


52.48 
36.49 


105.76 
103.92 


KLSQ 
KLSE 


62.76 
62.93 


51.95 
54.10 


109.13 
103.27 




e(U f )% 


e(V f )% 


e{W f )% 




e(U f )% 


e(V f )% 


e(W f )% 


KLSQ 
KLSE 


5.08 
3.78 


19.40 
15.92 


122.54 
86.68 


KLSQ 
KLSE 


4.18 
3.61 


15.87 
16.50 


96.63 
72.35 








SC (c) 


e(U)% 


e(V)% 


e(W)% 




SC (d) 


e(U)% 


e(V)% 


e(W)% 


KLSQ 
KLSE 


10.36 
10.42 


38.47 
38.59 


112.72 
108.99 


KLSQ 
KLSE 


10.21 
10.37 


37.75 
38.19 


109.98 
108.10 




e(U')% 


e(V')% 


e(W')% 




e(U')% 


e{V')% 


e(W')% 


KLSQ 
KLSE 


65.02 
65.37 


54.55 
54.70 


110.30 
106.68 


KLSQ 
KLSE 


64.15 
65.08 


53.54 
54.12 


107.56 
105.82 




e(U f )% 


e(V f )% 


e(W f )% 




e(U f )% 


e(V f )% 


e(W f )% 


KLSQ 
KLSE 


5.61 
4.22 


22.97 
18.51 


137.30 
97.95 


KLSQ 
KLSE 


5.02 
4.03 


20.16 
16.41 


121.39 
91.39 



SC (e) 


e(U)% 


e(V)% 


e(W)% 


KLSQ 


11.30 


40.99 


119.15 


KLSE 


10.24 


38.35 


107.59 




e(U')% 


e(V')% 


e(W')% 


KLSQ 


71.01 


58.17 


116.48 


KLSE 


64.23 


54.31 


105.33 




e(U f )% 


e(V f )% 


e(W f )% 


KLSQ 


6.78 


25.30 


166.46 


KLSE 


3.70 


16.61 


87.30 



Table 4: Relative percent errors, in L 2 norm, in the reconstruction of the velocity components (e(U), 
e(V), e{W)\ of their fluctuating part (e(U' ) , e(V ) , e(W' )), and of the part of the velocity field 
projected on the retained POD modes (e(U /), e(Vf), e(Wf)). The errors are averaged in time in the 
estimation interval, which starts just after the calibration interval. SC: Sensor Configuration. 



The parameter Cr of equation [8] was selected by experimenting different values. For example, 
in figure |4] we show the L 2 relative error in the reconstruction of the U component projected on the 
retained POD modes, as a function of Cr. The results are relative to the configuration (b) in the time 
interval close to the calibration one, using K-LSE. For all of the results shown in the following we 
took Cr = 0.1 for K-LSE. Note that for Cr > 10 2 the results are basically those of a simple LSE. 
A similar analysis was performed using K-LSQ and the optimal value that was selected is Cr = 10. 



INRIA 



A non-linear observer for unsteady three-dimensional flows 



17 



Results relative to the configuration (b) in the time interval close to the calibration one are reported 
in figure [5] where some representative modal coefficients predicted by the calibrated POD model, 
LSQ, LSE, K-LSQ, K-LSE and SLSE are plotted, together with the projection of the DNS velocity 
fields on the corresponding POD mode. 

Results of configuration (b) has been shown because the placement of the sensors is appropri- 
ate for the LSE method, as already discussed, and this makes the comparison with the proposed 
approaches more comprehensive. The results obtained with the other sensor configurations and on 
the time interval far from the calibration one are qualitatively analogous and quantitatively similar 
to those reported in figure [5] except for the LSE and the LSQ methods which are more sensitive 
to sensor placement. In figure [5] it is seen that LSE, LSQ, and also the calibrated POD model 
provide reasonable predictions only for the first two modal coefficients, that are associated with 
the vortex-shedding dynamics. The predictions of the remaining modes are completely unreliable. 
When dynamic estimation is applied, or when the SLSE approach is used, predictions are definitely 
improved. In particular, this is true for modes like oi or an that are related to the vortex-shedding, 
i.e., almost periodic with a period that is the same or a multiple of the vortex-shedding period. The 
prediction of the remaining modes is definitely less accurate (see 03 and 020), especially when very 
low frequencies are dominant, as in the case of a 3 . However, the overall accuracy is significantly 
improved in comparison with the LSQ and LSE approaches alone. The K-LSE approach is sys- 
tematically more accurate than SLSE. This is quantitatively confirmed by the relative error between 
the DNS velocity components and those reconstructed by the dynamic approaches and the SLSE 
one, reported in table [4] for all the sensor configurations. However, the table shows that there are 
no significant differences in accuracy between the K-LSQ, K-LSE and SLSE methods, and the re- 
construction errors can be considered satisfactory if we consider the complexity of the flow with 
respect to the 2D case. Indeed, the streamwise and vertical velocity components are reconstructed 
with errors of the order of 9% and 33%, respectively. Note that errors are computed over the whole 
computational domain, which extends for 20 L behind the cylinder, and, as it is qualitatively shown 
in figure|6] the reconstruction errors are small in the near-wake and they progressively increase mov- 
ing away from the cylinder in the downstream direction. Concerning the spanwise component of the 
velocity, errors are large since the retained POD modes themselves poorly represent this component 
of velocity as it is not energetically significant, in average, with respect to the remaining ones. This 
aspect might be improved working on the construction of the POD basis choosing, for instance, a 
different norm which weights more the spanwise component of the velocity or which corresponds to 
a quantity different from kinetic energy. The K-LSQ, K-LSE and SLSE methods are also similar in 
the sensitivity of the predictions to sensors type and placement, which is generally low. Neverthe- 
less, the predictions given by the K-LSE method are systematically the most insensitive to sensors 
placement. Table [4] also shows the relative reconstruction errors on the fluctuating part of the ve- 
locity components. The errors in the reconstruction of the fluctuating part of both the streamwise 
and the lateral velocity components are comparable. To help the interpretation of the errors in the 
reconstruction of the fluctuating flow field, which might seem large at first sight, it is important to 
point out that, outside the calibration interval, the amount of fluctuating energy that can be recovered 
by 20 modes is slightly more than 60% (0). Therefore, the error in the velocity components is due 
in part to a pure approximation error. In other words, the accuracy of the best possible reconstruction 
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SC (a) 


e(U)% 


e(V)% 


e(W)% 


KLSQ 
KLSE 


8.79 
8.96 


32.56 
33.54 


114.51 
110.95 




e(U')% 


e(V')% 


e(W')% 


KLSQ 
KLSE 


55.91 
56.97 


46.13 
47.53 


111.51 
108.00 




e(U f )% 


e(V f )% 


e{W f )% 


KLSQ 
KLSE 


4.13 
3.63 


15.74 
14.90 


101.87 
84.19 




SC (c) 


e(U)% 


e(V)% 


e(W)% 


KLSQ 
KLSE 


9.47 
9.17 


33.88 
33.19 


120.45 
110.95 




e(U')% 


e(V')% 


e(W')% 


KLSQ 
KLSE 


60.08 
58.36 


47.97 
47.03 


117.20 
108.09 




e(U f )% 


e(V f )% 


e(W f )% 


KLSQ 
KLSE 


5.28 
3.07 


19.02 
13.56 


133.11 
74.63 



SC (b) 


e(U)% 


e(V)% 


e(W)% 


KLSQ 
KLSE 


9.50 
8.83 


34.47 
33.11 


121.52 
111.35 




e(U')% 


e(V')% 


e(W')% 


KLSQ 
KLSE 


60.50 
56.17 


48.87 
46.91 


118.35 
108.42 




e(U f )% 


e(V f )% 


e(W f )% 


KLSQ 
KLSE 


5.15 
3.34 


19.82 
14.68 


118.19 
80.84 




SC (d) 


e(U)% 


e(V)% 


e(W)% 


KLSQ 
KLSE 


9.42 
9.44 


35.45 
34.40 


121.86 
110.40 




e(U')% 


e(V')% 


e(W')% 


KLSQ 
KLSE 


59.93 
59.93 


50.26 
48.71 


118.64 
107.54 




e(U f )% 


e(V f )% 


e(W f )% 


KLSQ 
KLSE 


5.11 
3.62 


20.46 
15.74 


133.56 
76.05 



SC (e) 


e(U)% 


e(V)% 


e(W)% 


KLSQ 


10.60 


38.50 


136.53 


KLSE 


9.22 


33.49 


110.24 




e(U')% 


e{V')% 


e(W')% 


KLSQ 


67.42 


54.59 


132.86 


KLSE 


58.65 


47.45 


107.38 




e(U f )% 


e(V f )% 


e(W f )% 


KLSQ 


7.29 


26.00 


183.25 


KLSE 


3.27 


14.10 


77.47 



Table 5: Relative percent errors, in L 2 norm, in the reconstruction of the velocity components (e(U), 
e(V), e{W)\ of their fluctuating part (e(U' ) , e(V ) , e(W' )), and of the part of the velocity field 
projected on the retained POD modes (e(U /), e(V/), e(Wf)). The errors are averaged in time in the 
estimation interval, which is about 52 time units far from the calibration interval. 



is limited from above by the capability of the POD modes to actually represent the flow outside the 
time interval where the snapshots were taken, which however, increases using a larger snapshots 
database, as shown in 0. 

In order to separate the accuracy of the estimation methods from the representativeness of the 
POD basis, we evaluated the errors between the estimated field and the projections of the velocity 
fields on the retained POD modes. Such errors are reported in tables [4] and [3] These errors are 
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(c) 

Figure 6: Isosurfaces of the velocity components u (left, grey = 0.5, dark grey = 1.0), v (center, grey 
= -0.25, dark grey = 0.25) and w (right, grey = -0.075, dark grey = 0.075) of a snapshot outside the 
database: (a) actual snapshot, (b) snapshot projected on the retained POD modes, (c) reconstructed 
snapshot using the K-LSE technique with the sensor configuration (b). 

systematically and significantly lower than those observed in the reconstruction of the whole flow 
fields. 

Comparing the results reported in table [4] with those of table [5j where the same quantities are 
shown for the time interval far from the calibration one, it can be seen that the accuracy of the 
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predictions is independent of the distance from the calibration interval, this being a very positive 
feature of all the proposed approaches. 

In figure [6] the velocity components obtained by DNS at t = 426.6 (a snapshot outside the 
database used for the derivation and calibration of the POD model) are plotted together with their 
projection in the space of the retained POD basis, which represents the best approximation of the 
flow which can be estimated with the retained POD modes, and with the prediction given by the 
K-LSE method. It can be seen that the main structures characterizing the streamwise and lateral 
velocity fields are well reconstructed. As for the spanwise velocity component, the reconstruction 
accuracy is not satisfactory, but this is due to the fact that it is one order of magnitude lower than the 
other components, as already discussed. 

5 Conclusions 

We devised a method to construct a non-linear observer for unsteady flows. This method is based on 
the coupling of a non-linear low-dimensional model of the flow with a linear technique that estimates 
the coefficients of the flow representation in terms of POD modes. The underlying idea is that the 
estimated flow should approximately satisfy the POD model. The coupling leads to a nonlinear 
minimization problem solved by a pseudo-spectral approach and a Newton method. 

The non-linear observer was applied to the laminar flow around a confined square cylinder at two 
different Reynolds numbers; at the first the flow is two-dimensional, while in the second case three- 
dimensional phenomena occur in the wake. In the two-dimensional case, since the flow patterns 
are rather simple, the results show that the proposed procedure is able to give a significantly more 
accurate estimate of the POD coefficients even with a limited number of sensors, than those obtained 
with the LSQ and the LSE approaches. The QSE method improves the predictions of the LSE, even 
if the overall accuracy remains definitely lower than those given by the proposed approaches. In 
the three-dimensional case, the flow dynamics is more complex, and not only LSE and LSQ, but 
also the calibrated POD dynamical system provide unsatisfactory coefficient estimations when used 
outside the calibration interval. Conversely, the proposed procedure, combined with either LSQ or 
LSE, gives accurate predictions of the coefficients of those POD modes that are related to vortex 
shedding. For the remaining modes, the accuracy is lower. Nevertheless, the instantaneous velocity 
fields are reconstructed with satisfactory accuracy, both close and far from the calibration interval 
of the POD model. Moreover, K-LSE and K-LSQ methods are weakly sensitive to sensor type and 
placement. The results obtained with the proposed approaches are comparable to those obtained 
by the SLSE approach, which also uses the temporal history of the flow measurements, but in the 
Fourier space. However, this latter techniques has computational complexity which is significantly 
larger than that of LSE or LSQ, and it is comparable to that of the proposed approaches. 

This work was funded in part by the HPC-EUROPA project (RII3-CT-2003-506079). IDRIS 
(Orsay, France) and M3PEC (Universite Bordeaux 1, France) provided the computational resources. 
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